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SUMMARY 


Slender-body aerodynamic theory is used in conjunction with thin- 
plate theory in the flutter analysis of low-aspect-ratio rectangular 
wings of constant thickness when chordwise variations of deflections 
are considered. The spanwise variation of deflection is given by a 
parabola, and the chordwise variation is allowed complete freedom. The 
results show the variation of flutter speed and mode shape with aspect 
ratio. Comparisons are made with results obtained by approximating the 
chordwise deflection shape by the first few terms of a power series . 
Comparisons with some preliminary experimental results are also included. 


INTRODUCTION 


The prediction of flutter of wing and tail surfaces of very low 
aspect ratio is a problem of some concern to aircraft designers . The 
difficulties in the flutter analysis of such surfaces are mainly con- 
nected with the presence of large amounts of chordwise curvature in the 
flutter mode . It is of interest to investigate the complexity of the 
chordwise deflection shape at flutter and to determine to what degree 
of accuracy the chordwise deflection must be represented in order to 
obtain good results . This paper is concerned with the flutter behavior, 
both theoretical and experimental, of the simple low-aspect-ratio con- 
figuration shown in figure 1. The analysis is similar t.o that of ref- 
erence 1, which treated the static divergence behavior of the same 
configuration . 


SYMBOLS 


x,y,z coordinate system (see fig. l) 

w(x,y,T) wing deflection, positive in z-airection 
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free -stream velocity 
wing thickness 
wing chord 
wing semi span 
time 

chordwise deflection shape 


dynamic pressure. 


PV £ 


Young's modulus of elasticity 
free-stream density of fluid 
density of material 
period of oscillation 
Mach number 

perturhation-velocity potential 
flutter-speed parameter, ^(l - ^ -S— 


Pmts^ 


flutter-frequency parameter, ca 

mass-ratio parameter, £ 

^ Pm * 3 

Poisson's ratio (taken as 1/3 in all computations) 

plate stiffness in bending, Et^/ 12 (l - p.^) 
flutter frequency 


THEORETICAL APPROACH 


The configuration analyzed in this paper consists of a rectangular 
plate of constant thickness cantilevered from a rigid wall. (See fig. 1.) 
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This plate may be thought of as representing one -half of a wing with a 
chord c and a semispan s . The plate is located in a fluid flow with 
a free- stream velocity V. The deflection shape of such a low-aspect- 
ratio plate can be expected to vary in a much more complicated manner 
in the chordwise direction than in the spanwise direction- For this 

reason the deflection w is assumed to vary as w(x,y,T) = y^F(x,*r), 
where the spanwise deflection is given by a simple parabola and the 
chordwise variation F of the deflection is an arbitrary function of 
the chordwise coordinate x and time t. The distortions of the plate 
are found through the use of ordinary thin-plate theory. The aerodynamic 
loadings are found most simply by using slender-body aerodynamic theory. 
In this theory streamwise perturbations are neglected in comparison with 
perturbations in the crossflow direction. The use of this approximate 
aerodynamic theory simplifies the aeroelastic problem to the extent that 
an exact solution is possible. A brief description of the analysis and 
the resulting equations are given in the appendix. 


RESULTS AMD DISCUSSION 

Some results are given by the boundaries shown in figure 2. The 
ordinate is the dynamic-pressure parameter SJB. — } in which q. is the 

(t/s) 5 

dynamic pressure, E is Young* s modulus of elasticity for the plate, 
and t/s is the ratio of thickness to semispan. The abscissa is the 
ratio of chord to semispan c/s. The variation of the dynamic-pressure 
parameter for flutter with the ratio of chord to semispan is dependent upon 

the mass -ratio parameter in which -2- is the ratio of air density 

Pm"k Pm 

to plate density. The flutter boundaries are given for two values of 
mass-ratio parameter; the region above a particular boundary is unstable 
while that below is stable. Also shown is the result for static diver- 
gence obtained from reference 1. This result, which is, of course, 
independent of the mass ratio, is indicated by a single curve. Note 
that divergence is less critical than flutter for these particular mass 
ratios. Note also that, for the higher mass ratio, the flutter boundary 
consists of a series of loops approaching a constant value of dynamic- 
pressure parameter. The lower curve also has this characteristic; how- 
ever, the loops are so elongated in this case that only one can be seen 
in this figure. 

The kind of flutter mode shapes obtained from the analysis is 
shown in figure 3 . The top set of curves shows the components of tip 
deflection which are in phase and out of phase with the maximum leading- 
edge deflection for a value of c/s corresponding to the tick mark on 
the bottom of the first loop on the flutter boundary in figure 2. The 
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tip deflection is given when the leading edge has its maximum amplitude 
and at one-quarter of a period later when the leading edge has zero 
deflection. The bottom set of curves gives the components of mode shape 
for a value of c/s given by the tick mark on the second loop of the 
flutter boundary of figure 2. The effect of increasing the chord is to 
add more waves to the mode shape . 

In the results discussed so far, the chordwise variations of deflec- 
tions were allowed complete freedom and an exact solution was possible. 

In a practical case, an exact solution would not be feasible, and some 
sort of approximation of the chordwise deflection shape would be nec- 
essary. Some flutter boundaries obtained by approximating the chord- 
wise deflections by the first few terms of a power series are shown in 
figure 4. The dashed curve gives the results for parabolic deformations, 
and the long-and- short-dashed curve gives the results for cubic deforma- 
tions. The exact boundary is also shown for comparison. Both approxi- 
mations yield good results for the lower values of c/s. The cubic 
approximation is almost exact. For longer chords, however , both approxi- 
mations yield poor results . Apparently, in order to analyze the flutter 
behavior of wings in this range, higher order terms in the deflection 
shape must be used. 

It can be seen from figure 2 that for each value of mass ratio a 
limiting value of the dynamic-pressure parameter can be obtained by con- 
sidering the values at the bottom of the loops as the chord becomes 
large . The variation of this limiting value with mass ratio is shown 
in figure 5 • The results obtained so far indicate that this curve gives 
a conservative estimate of the flutter speed for thin rectangular plates 
of very low aspect ratio. It should be noted that the flutter speed is 
less than the divergence speed and seems to approach it asymptotically 
for high mass ratio. 

Several preliminary tests of some low-aspect-ratio plates at super- 
sonic Mach numbers have been run in the Langley 9- "by 18-inch supersonic 
flutter tunnel. The results of these tests are shown in figure 6. The 
flutter speed of the various models is shown as a ratio of experimental 
flutter speed to calculated flutter speed, where the calculated flutter 
speed was obtained from the curve in figure 5* This ratio is plotted 
against the ratio of chord to semispan. The agreement between theory and 
experiment is fairly good in view of the approximations inherent in the 
theory and of the preliminary nature of the tests . 

During the tests the flutter modes were observed, and it was noted 
that the specimens with the larger ratios of chord to semispan had the 
more complicated mode shapes as is predicted by theory. 
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CONCLUDING REMARKS 


The flutter analysis of low-aspect-ratio rectangular plates indi- 
cates that the flutter-mode shape has an increasing number of waves in 
the chordwise direction as the aspect ratio is reduced; approximating 
the chordwise deflection shape by parabolic or cubic curves yields flut- 
ter speeds in fair agreement with those of the more exact theory, pro- 
vided that the aspect ratio is not too low. The cubic curve gives some- 
what better results . For lower aspect ratios, higher order approximations 
must be used. Experimental data indicate that the flutter speed and the 
type of mode shape yielded by the theory are in fair agreement with 
experimental results . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 7 , 1957 - 
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APPENDIX 

ANALYSIS 


In the application of slender -body aerodynamic theory to the present 
unsteady-flow problem, terms containing time derivatives in the velocity- 
potential equation for linearized flow are neglected in addition to those 
containing streamwise x-derivatives . The velocity-potential equation 
thus reduces to Laplace's equation in the crossflow plane. The boundary 
conditions on velocity and the pressure-potential relations are taken to 
be the same as those ordinarily used in unsteady linearized aerodynamic 
theory . 

The potential o and the aerodynamic loads resulting from the 
given deformation shape are calculated in a manner similar to that used 
in reference 1 for the static-divergence problem. 

The principle of minimum potential energy is used to derive the 
differential equation of equilibrium for the function F(x, t) in a man- 
ner analogous to that used in reference 1. Solution yields an eigen- 
value equation which relates the flutter speed to the properties of the 
plate and the surrounding air: 


0 = 64<x3 cos + Et_ cosh Sjp cosh ^ sinh S& sinh + 


i(lu sin cosh &=• + cosh sinh 

B Ok S Sp S S / 


( 1 ) 


where 
i = \[-l 


D 1 = “[^1^ + A 1 B 2 + A #lJ 

E 1 = ( A 3 + >B 5 )^2A 5 - (by 2 ' + a 2 + p 2 ) Aj^j + 
(A4 - >,)[a*, 2 p 2 A4 - (4/ 2 + a 2 + 3 2 }aJ 
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F x = (Aj + 7 vB 5 )^ 2 a 2 p 2 A 4 - (V 2 + a 2 + P 2 )aJ + 
a 2 p 2 (A^ - A)j^ 2 A 5 - (A 7 2 + a 2 + P 2 )a^J 
d 2 = “( A 1 B 2 " A 2 B l) 

xfsr 2 + Sa-L - a\ + ff ) - 2(r 2 - A ) £ 
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A2 = 


^ (1 - |i)(l - 2 H) - ^ (1 - jx)A + ^ (1 - H) 7 2 + 
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\ = & (1 - n) + 7 2 



In equation (l), a and |3 can be either real or imaginary quantities. 
Solution hy trial yields the curves given in figure 2. 

The divergence "boundary obtained from reference 1 can also be 
obtained from the present analysis by setting the flutter -frequency 
parameter K equal to zero. 


Mode Shapes 

The flutter -mode shapes shown in figure 3 can be obtained from the 
following equation, where the real part gives the component in phase and 
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the imaginary part, the component out of phase with the maximum leading- 
edge deflection: 


f(x) = 


a-e 1 * 1 s + Cge^ 2 s + C*e 




+ Cji.e 


Cn + Co + C' 


3 + ci, 


where 


Cl = (e. 2 - m^jA^e 1 ( m l' Ha 3 ) s 4 

(*, - ■*,)(*, - + 

)(a 2 -^ e -H»l^)w] 


l m 4 • m 3 


C 2 


= [(m 3 - mi,) (a 2 - > 1 + 

( m 4 - m l) ( a 2 - ®3 2 ) A i4 e 1 ( m2+in ' 5 ^ 8 + 
( m l " m 3) (*2 “ m 4 2 ) A l3 e ^ S J 


C 3 “ 


(m 4 - mg) (a 2 - m^A^e^l 4 ^) I + 
( mi - mi,) ^a 2 - m 2 2 )A 1 i t e 5 ^ 8 + 

(mg - m 1 )(a 2 - mi, 2 )A 12 e 1 ^ 3+ ^ ® ] 
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c 4 = [(“2 - m 3)( a 2 - “ l 2 )^ 6 ^ 1+4 ^ s - 

( m l - “2)(a2 - m 3 2 ) A 12 e 1 ( m 3 +m 4 ) .] 


in which 


m-L = 7 + ip 

»2 = 7 - 10 
m j = -7 + la 
nj, = -7 - la 

A pq = -^a 2 a 5 - ag^, 2 + m-pin^ + m q 2 ) - 4a 5 iynq + nip 2 !]^ 2 


*2 = y M- 

a = 10 J 

a 3 3 \ 4 *7 


Approximate Solutions 

The approximate analysis used parallels the analysis for divergence 
presented in reference 1. Cubic approximation of the chordwise varia- 
tion of deflection results in the following complex determinant: 



I 


[* < i2K \|^] [$ A + VK f + 

M [j ik ^ (i - >*) ?+ 

Ei * " “ 11 Cl [i A + T( ,t -5">5" i f K 

[i A + ion fj l A + ^(2-(i)f ,.1 kJ§ 


where 


^ a + (o\ + |i n) g i. 

^ A + M (4 - jn)J S ,• 1 fl K^| 

jA-jf (*-»)}* + 4^ + iacj^ 

^A + |a + ^(3-2,ojf + S^ + ifiKj§ 


A = 


20 - K®(l + g 


[jj A + (12X + 10n) f |. illK^S 
1 A + |o\ + as (2 - |»)j fi + i?K ^ 
5 a + j& + a (i - 2n)| f * 6 a| + i is X 

[7 A H ( i f * * 0(5 ■ 2|1) ] f ■"• M ij + 12K 


' 0 (2) 


r Ihis equation was solved by trial to obtain the long-and-short-dashed curve in figure 4 . The 
equation for parabolic deformation is obtained from equation (2) by omitting the last column 
and last row. The results for parabolic deformation are shown by the dashed curve of figure 4 . 
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FLUTTER MODE SHAPES 



Figure 3 


FLUTTER BOUNDARIES OBTAINED FROM APPROXIMATE 
DEFLECTION FUNCTIONS 



Figure 4 







